Nonadiabatic charge pumping in a one-dimensional system of noninteracting electrons 

by an oscillating potential 
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Using a tight-binding model, we study one-parameter charge pumping in a one-dimensional system 
of non-interacting electrons. An oscillating potential is applied at one site while a static potential is 
applied in a different region. Using Floquet scattering theory, we calculate the current up to second 
order in the oscillation amplitude and exactly in the oscillation frequency. For low frequency, the 
charge pumped per cycle is proportional to the frequency and therefore vanishes in the adiabatic 
limit. If the static potential has a bound state, we find that such a state has a significant effect on 
the pumped charge if the oscillating potential can excite the bound state into the continuum states 
or vice versa. Finally, we use the equation of motion for the density matrix to numerically compute 
the pumped current for any value of the amplitude and frequency. The numerical results confirm 
the unusual effect of a bound state. 

PACS numbers: 73.23.-b, 73.63.Nm, 72.10.Bg 
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I. INTRODUCTION 

The phenomenon of charge pumping and rectification 
by time-dependent potentials applied to certain points in 
a system has been extensively studied both theoretically 
[1-31] and experimentally [I| Q M, M, M, El . The 
idea of charge pumping is that periodically oscillating 
potentials can transfer a net charge per cycle between two 
leads which are at the same chemical potential. For the 
case of non- interacting electrons, theoretical studies of 
charge pumping have used adiabatic scattering theory 0, 
Hoi . IrH [l2| , Floquet scattering theory 0, [l7| , variations 
of the non-equilibrium Green function formalism [l8l . Il9l . 
|20| , and the equation of motion approach [3l[ . The case 
of interacting electrons has also been studied, using a 
renormalization group method for weak interactions [391 ] , 
and the method of bosonization for arbitrary interactions 
[40-50]. 

Apart from a few papers [H, H(| [H[ , the earlier stud- 
ies of charge pumping have generally studied systems in 
which oscillating potentials are applied to two or more 
sites. In such cases, it is known that if the oscillation 
frequency lo is small, the dc part of the pumped current 
is proportional to to; the charge pumped per cycle (with 
time period 2tt/lo) therefore has a finite value in the adi- 
abatic limit lo — > 0. However, an oscillating potential 
applied to a single site can also pump charge provided 
that the system has no left-right symmetry, as has been 
emphasized in Refs. [l9l [20| . This can happen if, for 
instance, there is a static potential at another site. In 
these cases, however, the dc part of the pumped current 
is proportional to lo 2 , and the charge pumped per cycle 
is proportional to a; if a; is small. In this paper, we will 
study such a system in detail using both analytical and 
numerical methods. Our analysis will not be restricted 
to small values of lo. If lo is large enough, we discover 
that a bound state (defined as a state whose energy lies 
outside the continuum band of the tight-binding model 
that we will consider) can have a significant effect on 



the pumped charge. Namely, although the bound state 
has a localized wave function and therefore cannot con- 
tribute to the current at the same energy, the electrons 
can be scattered from the bound state to the continuum 
states (lying within the band) or vice versa, and this can 
contribute to the current flow. This phenomenon does 
not seem to have received much attention in the existing 
literature. 

The plan of the paper is as follows. In Sec. II, we will 
use Floquet scattering theory to obtain an expression 
for the pumped current produced by a single harmoni- 
cally oscillating potential, when there is a static potential 
present at some other point in the system. Our analysis 
will be exact in the scattering matrix arising from the 
static potential and in the oscillation frequency, but it 
will be pcrturbative in the amplitude of the oscillations. 
The effect of a bound state produced by the static po- 
tential will also be considered using the same formalism. 
In Sec. Ill, we will use the equation of motion to numer- 
ically compute the pumped charge for the same model; 
this method can be used for any value of the oscillation 
amplitude. Our numerical results will confirm the un- 
usual effect that a bound state can have on the pumped 
charge, as well as the difficulty which the equation of mo- 
tion method faces in dealing with a bound state [5l[ . We 
will summarize our results in Sec. IV. 



II. FLOQUET SCATTERING THEORY 



Let us consider a one-dimensional system consisting of 
two semi-infinite leads a — L, R (denoting left and right) 
and a finite region called the wire which lies between the 
two. We will model the three regions together by a lat- 
tice model of spinless electrons governed by a one-channel 
tight-binding Hamiltonian with the same hopping ampli- 
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tude —7 on all bonds, namely, 

00 

Hq = -7 E ( cl +1 c n + c+Cn+i) + ac\c , (1) 

n— — 00 

where we have added a time-independent potential with 
strength a at the site n = 0. The dispersion of the elec- 
trons in the leads is Ek = —27 cos A:, where k lies in the 
range [— 7r, tt\. (We are setting the Planck constant H 
and the lattice spacing equal to unity). The two leads 
are assumed to be at the same chemical potential /x and 
temperature T. 

The static potential at the site n = causes scattering 
of electrons incident from the left or right lead. If an 
incident electron has wave number k, the effect of the 
potential is described by a scattering matrix 



S(k) 

where r R (k) 
t R {k) 



tR r R 



r L {k) 
t L (k) 



27 sin k + ia 
27 sin k 



27sinfc + ia 



(2) 



and fL(R) and ^.r(l) denote the reflection and transmis- 
sion amplitudes for an electron coming from the left 
(right) respectively. 

It turns out that for any value of the parameter a =/= 
0, there is a bound state. For a > 0, the bound state 
energy Eb = \f 47 2 + a 2 lies above the continuum, and 
the normalized bound state wave function is 



ipB(n) = (-1)" Vtanhrc e" KW for all n, (3) 

where n = sinh~ 1 (a/27). For a < 0, the bound state 
energy Eb — — y ij 2 + a 2 lies below the continuum, and 
the corresponding wave function is 



ipB(n) = V tanh k e 



for all n, 



(4) 



where k = sinh _1 (— a/27). 

We now apply an oscillating potential at the site n = r 
of the model described in Eq. ((T|), where we assume that 
r > 1. This part of the Hamiltonian is given by 



V(t) = b cos(wi) c\.c r . 



(5) 



Then Floquet scattering theory works as follows [16|, ll7| . 
Incoming electrons of energy Eq gain or lose energy in 
quanta of to on interacting with the oscillating potential. 
Hence, the outgoing states are characterized by energies 
E p = Eq+puj, where p = 0, ±1, ±2, ■ • • ; the energies with 
p / are called the Floquet side bands. The effect of the 
oscillating potential can be described by a Floquet scat- 
tering matrix S a p(E p , Eq), which is the amplitude for an 
electron with energy Eq entering through lead (3 to leave 
with energy E p through lead a. In the leads, the prop- 
agating modes have energies lying within the continuum 
band [—27, 27]; only these modes can directly contribute 



to charge pumping. States with energies lying outside 
the continuum band have wave functions which decay 
exponentially into the leads and hence do not directly 
contribute to charge transfer. The wave function of an 
electron coming from the left lead with an energy Eq and 
wave number fco (with Eq = — 27cosfco) is given by 

^(n) = e l(kon - Eot) + r ^,v e' ( -^ n - E ^\ (6) 

p 

at a site n far to the left of the scattering region, and 

tl>(n) = e l( ^"-^ 4) , (7) 

p 

far to the right of the scattering region, where E p = 
— 27 cos k p , and the sums over p run over values such 
that E p lies within the continuum band of the leads. The 
quantities rL,p and tji yP denote reflection and transmis- 
sion amplitudes in the different side bands; they respec- 
tively denote the elements Sll{E p , Eq) and Srl(E p , Eq) 
of the Floquet scattering matrix, where L and R denote 
the left and right leads. Similarly, the wave function of 
an electron coming from the right lead with an energy 
Eq and wave number fcp is given by 



^(n) = e i<.-*on-&>t) + 



p 



i(k p n—E p t) 



> (8) 



far to the right of the scattering region, and 



far to the left of the scattering region. Due to unitarity, 
we have the relations 



E - 1 KpI 



\tR,p\ 



1, 



and V * [ \r R .p\ 2 + \t L J 2 ] = 1, (10) 



where v p = 27 sin k p is the velocity in the p-th side band. 

The reflection and transmission amplitudes are found 
by writing down the wave functions in the scattering re- 
gion, and matching coefficients of terms having the same 
time dependence (e ±tEpt ) in the Schrodinger equation at 
different sites. If the oscillating potentials are weak, the 
reflection and transmission amplitudes decrease rapidly 
as \p\ increases; at first order in the potentials, only 
p = ±1 contribute. The dc part of the current in, say, 
the right lead is then given by 



cLEq Vl 2 



FR,i 



x {f(E ) - f(E 1 )} 



27 



-2 7+W 27r Vq 

x {f(E ) -/(£?_!> ], (11) 
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where f(E) = \j\^ E ~^l kBT + 1] is the Fermi function, 
and q is the charge of an electron. The upper limit is 27— 
w in the first integral in Eq. (fTTj) because ii E > 27 — oj, 
Ei will lie above the continuum band and will therefore 
not contribute to the current. Similarly, the lower limit 
is — 27 + oj in the second integral in Eq. (fTTj) because if 
Eq < —2"f + oj, E-i will lie below the continuum band 
and will not contribute to the current. If oj > 47, the 
integrals will contribute nothing since all the side bands 
lie outside the continuum band; the pumped current will 
therefore vanish. 

For the case of a single static potential and a single 
oscillating potential described by Eqs. ([p} and (O, we 
find that 



ib 



l R,l 



47 sin fci 
x [ 1 + r fl (fei) e l2kir 



e - l (fco+fc 1 )r j j + e i2kor 



ib 



t R (k ) e l{k ^ r 



47 sin ki 
x [ 1 + Pfl (fci) e i2kir 



(12) 



where the functions r R and tu are given in Eq. ([2]) . The 
expressions for rfj._i and t R ^-\ can be obtained from 
Eq. (fT2"|) by replacing k\ by k-%. By substituting all 
this in Eq. (jlip . we obtain an expression for I R which 
is valid up to second order in the dimensionless quantity 
6/(27 sinfcp). 

In the limit oj — ► 0, we have f(E±i) — f{E a ) 
■ ojdf(E )/dE . By shifting E — ► E + oj in the sec- 
ond integral in Eq. (fTTj) . one can show that Ir goes 
as oj 2 for oj — ► 0. In particular, for zero temperature, 
df (Eq) I dE = —5(E — fi), and we get 



qoj 2 b 2 \t R (k F ) 



d 

7k- 



|1 + r R (k) e 



i1kr\2 



647r 7 3 sin 3 fc F \dk \t R (k)\ 2 , k=kf 

(13) 

where kp is the Fermi wave number given by jjb = 
— 2 j cos k F - Note that Eq. (fT3"|) vanishes if r = since 
1 + r R (k) = t R (k); thus a non-zero value of r, i.e., a 
left-right spatial asymmetry, is necessary to have charge 
pumping. The expression in (| 1 3[) is in contrast to the 
case in which there are two or more oscillating potentials 
at different locations in which case the pumped current 
generally goes as oj for oj — > 0. 

Equation (fTB")) can be generalized to the case where 
the time-dependent potential at a given site has a finite 
number of oscillating terms each of which has a small 
amplitude and low frequency, namely, if bcos(oji) in Eq. 
© is replaced by bi cos(ojit + To second order 
in bi and oji, we find that the dc part of the pumped 
current at zero temperature can be obtained by replacing 
the factor oj 2 b 2 in Eq. (J3]) by 

We have so far discussed the contribution to the 
pumped current from the scattering states only. It turns 
out that a bound state can also contribute to the pumped 
current. Suppose that the static potential gives rise to 
a bound state with energy E R < —2j- If the oscillation 



frequency is large enough that E\ — E R + oj lies within 
the continuum band [—27, 27], then the bound state will 
contribute to the current. Using Floquet theory up to 
first order in b, we find that the wave function is given 
by a R e l{kin ~ Elt ^ and a L e l( ~ kin ~ Elt ^ far to the right and 
left of the scattering region respectively, where 



ib 



a R 



— ik\r 



47 sin ki 
x [ 1 + r R (h) e t2k ' r }, 



ib 



47 sin k\ 



e lkir ^ B (r) t L {kx), (14) 



and the function is given in Eq. ([4]). The dc part of 
the current pumped to the right is then given by 

Irb - q 2 7 sinfc 1 [ \a R \ 2 - \a L \ 2 } {f(E B ) - f{E{)}, 

(15) 



where 27sinfci is the electron velocity. Equation (fT5|) 
has to be added to Eq. (fTTj) in order to obtain the total 
pumped current. At zero temperature, I RR is non-zero 
only if Eg + oj > [i > E R , and it is then independent of 
fj,. 

Similar considerations hold if there is a positive en- 
ergy bound state with E R > 27 and E-\ = Eb — oj lies 
within the continuum band. Such a bound state will then 
contribute to the pumped current. One can compute 
this contribution by applying a particle-hole transforma- 
tion to the calculation given above for a negative energy 
bound state. Under the transformation c n — > (— l) n c ra , 
the hopping term in Eq. ([1]) remains the same but 



c n c n 



changes sign. This is equivalent to changing a — > —a in 
Eq. ([T]), and the chemical potential /i — * — /i. Thus the 
filling fraction / = kp/it changes as / — > 1 — /, and the 
current changes as I R — > —Ir\ the latter can be seen di- 
rectly from the form of the current operator given in Eq. 
([IT]) below. 

We now observe that, to second order in the ampli- 
tude 6, Eqs. (|1H12|) and Eqs. (|14I15|) for the pumped 
current remain valid for the case of a general static po- 
tential which may be extended over more than one site, 
as long as the oscillating potential lies outside the region 
of the static potential; we only have to use the appropri- 
ate expressions for the reflection and transmission ampli- 
tudes r R ,±\(k) and tR t ±i(k), and the bound state wave 
function %1>b[t). The same derivation which was used to 
obtain the above equations for the case of a static po- 
tential at one site will work for a more general case. We 
can show this in a different way by using the idea of sum 
over paths. To first order in b, the reflection amplitude 
r R .\ in Eq. (|12p can be understood as arising from the 
sum of the following four paths (see Fig. 1), remember- 
ing that the oscillating potential lies to the right of the 
static potential. An electron incident from the far right 
with a wave number fco can be 

(i) reflected to the right by the oscillating potential to a 
wave number k\ with an amplitude —16/(47 sin ki) [this 
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a 



b 



ii) 



a 



b 



Hi) 



b 



iv) 



a 



b 



FIG. 1: (Color online) Different paths contributing to rn,i at 
first order in b; a and b denote the static and oscillating poten- 
tials respectively. The crosses denote the scattering events. 



amplitude can be derived using the Born approxima- 
tion on a lattice [52| and remembering that bcos(u>t) — 
(b/2)(e iuJt + e~ iuJ % or 

(ii) reflected to the right by the static potential with an 
amplitude r^(fco), and then transmitted to the right by 
the oscillating potential to a wave number k\ with an 
amplitude — i6/(47sinfci), or 

(iii) transmitted to the left by the oscillating potential 
to a wave number k\ with an amplitude — i6/(47sinfci), 
and then reflected to the right with an amplitude rn(ki), 
or 

(iv) reflected to the right by the static potential with an 
amplitude rn(ko), then reflected to the left by the oscil- 
lating potential to a wave number k± with an amplitude 
—16/(47 sin fci), and finally reflected to the right by the 
static potential with an amplitude rn(ki). 

Similarly, the transmission amplitude t^i in Eq. (112[) 
can be understood as a sum of the following two paths. 



An electron incident from the far left with a wave number 
ko can be 

(i) transmitted to the right by the static potential with 
an amplitude tn(ko), and then transmitted to the right 
by the oscillating potential to a wave number k\ with an 
amplitude — i6/(47sinfci), or 

(ii) transmitted to the right by the static potential with 
an amplitude tn(ko), then reflected to the left by the 
oscillating potential to a wave number ki with an ampli- 
tude —16/(47 sin fci), and finally reflected to the right by 
the static potential with an amplitude rn(k\). 

Similar ideas can be used to derive the expressions in 
Eq. (|14[) . The amplitude an is the sum of two terms. 
An electron with a wave function ipsir) can either be (i) 
transmitted to the right by the oscillating potential to 
a wave number ki with an amplitude — ib/ (Aj sin ki) , or 
(ii) transmitted to the left by the oscillating potential to 
a wave number k\ with an amplitude —16/(47 sin fci), and 
then reflected to the right by the static potential with an 
amplitude r^(ki). The amplitude ol corresponds to an 
electron with a wave function ipsif) being transmitted 
to the left by the oscillating potential to a wave number 
k\ with an amplitude — i6/(47sinfci), and then transmit- 
ted to the left by the static potential with an amplitude 

t L (h). 

We thus see that Eqs. (fl"2")l and (fl"4"|) are valid to first or- 
der in 6 for any static potential, provided that the static 
and oscillating potentials are separated by a finite dis- 
tance. 



III. EQUATION OF MOTION METHOD 

We will now discuss how the pumped current can be 
obtained by numerically solving the equation of motion 
for the density matrix of a system with a finite number 
of sites [3l|- The density matrix of the system evolves 
according to the equation of motion 



pit + dt) 



3 — iff(t)dt 



P(t) 



JH(t)dt 



(16) 



where H(t) = H + V(t) is given in Eqs. CQ) and ©. 
The current across any bond is then obtained by taking 
the trace of the appropriate current operator with p. The 
current operator on the bond from site n to site n + 1 
and its expectation value at time t are given by 

Jn+l/2 = iq"f ( c n+l c n ~ ci c n+l) , 

and j n+1 / 2 (t) = Tr{ p(t) j n+1 / 2 ) 

= i<n \p n<n+ i(t) - p n+ i. „(<)]. (17) 

The charge transferred between the left and right halves 
of the system L and R between two times can be found 
either by integrating the above expression in time, or by 
taking the operator 



AQ 



E 



J. 



E 



4, 



(18) 
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and computing Tr (p(t)AQ) at the two times; these 
methods give the same result for the charge transferred 
in a cycle. 

In all our calculations, we take the left and right leads 
to have Ni sites each and the wire in the middle to have 
N w sites; the total number of sites is N = 2Ni + N w . 
We take the density matrix at time t — to be given 
by that of a single system governed by the Hamiltonian 
Hq in Eq. ([1]) with N sites, chemical potential fi, and 
temperature T. If E a and ip a (n) are the eigenvalues and 
eigenstates of the Hq (a and n label the states and sites 
respectively), the initial density matrix is given by 

Pmn(0) = J2 ^(™) ^ ( 19 ) 
a 

We then evolve the density matrix in time and compute 
the current and charge transferred using Eqs. (|16til8|) . 

An important point to note is that the finite length 
of the leads (with TV; sites) implies that the system has 
a return time equal to 2Ni/vp where the Fermi velocity 
vp — 27sinfcp [5i(; this is the time required for an elec- 
tron to travel from the wire in the middle to the end of 
either of the two leads and then return to the wire. The 
numerical results can be trusted only for times which are 
less than the return time. Further, there are transient 
effects which last for one or two cycles; the effects of dif- 
ferent choices of the initial density matrix get washed 
out after this transient period. All the numerical results 
presented below are therefore taken from times which are 
larger than the transient time but smaller than the return 
time; typically, we have computed the charge transferred 
between the times 6tt/oj and lOir/u), where u> is the os- 
cillation frequency. The dc part of the charge pumped 
per cycle should of course be independent of the location 
of the bond where it is measured; we have checked that 
this is true for all our numerical results except for the 
contribution of a bound state as we will discuss below. 

We will now present our numerical results for the 
pumped charge (in units of q). In all cases we have set 
the temperature to zero, the hopping amplitude 7=1, 
and r — 5, i.e., the static and oscillating potentials are 
separated by five lattice spacings. 

In Fig. 2, we show the charge pumped per cycle (AQ = 
(27t/w)/r) as a function of the Fermi wave number kp for 
the case a = 1.5, b = 0.25, and uj = 7r/10 (corresponding 
to a time period of 20). The dash-dot and starred lines 
show the numerical and analytical results obtained from 
Eq. (fTT|) respectively; the agreement between the two 
is excellent. The pumped charge can be seen to go to 
zero at the band edges as expected. Although there is 
a positive energy bound state at Eb = 2.5, it does not 
contribute to the pumped charge because the first side 
band lies at an energy of £Li = Eb — u> — 2.19 which is 
above the continuum band. 

Figure 3 shows the charge pumped per cycle as a func- 
tion of kp for a = 1, b = 0.25, and uj = 7r/10. The dash- 
dot line shows the numerical results, while the starred 
and dash lines show the analytical results obtained from 
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FIG. 2: (Color online) Charge pumped per cycle from left 
to right versus the Fermi wave number for a static potential 
separated from an oscillating potential by five sites, with a = 
1.5, b — 0.25, and uj = n/10; the system has 252 sites. The 
numerical and analytical results are shown by dash-dot (red) 
and starred (blue) lines, respectively. 



Eq. (fTTj) (continuum states) and Eq. (fTTj) plus (fT5)) (con- 
tinuum and bound states) respectively. In contrast to the 
case shown in Fig. 1, the total pumped current does not 
go to zero near kp — it. This is because of the contribu- 
tion from a positive energy bound state; this has the en- 
ergy Eb — 2.24, and the first side band lies at an energy 
of = Eb — oj — 1-92 which lies within the contin- 
uum band. Hence, when the chemical potential exceeds 
E—i, the bound state begins to contribute to the pumped 
charge. The above value of E_\ corresponds to a Fermi 
wave vector of kp = cos _1 (— 1.92/2) = 2.86; we can see 
from the figure that when kp/ir exceeds 2.86/7T ~ 0.91, 
the total pumped charge begins to deviate from the con- 
tinuum contribution (which, according to Eq. lfTT|l , does 
go to zero as kp — > ir). Using Eq. (fT5]) . we can compute 
the contribution to the pumped charge from the bound 
state, AQ B = {2tt/uj)I rb . We find that AQ B = -0.0055 
which agrees reasonably well with the value obtained 
numerically when kp — > tt. However, we find that the 
pumped charge arising from the bound state has rather 
long lived transients, and the value of the dc part of 
the current varies significantly depending on the location 
of the bond where it is measured; the numerical result 
shown in Fig. 3 is the current measured at the bond ly- 
ing midway between the static and oscillating potentials. 
Thus our model has a numerical difficulty in correctly 
computing the contribution of a bound state to the cur- 
rent. The reason for this difficulty will be discussed in 
Sec. IV. 

Figure 4 shows the charge pumped per cycle as a func- 
tion of the oscillation frequency uj for a = 1.5, b = 0.25, 
and kp = 7r/4 (corresponding to fi ~ —1.414). The dash- 
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FIG. 3: (Color online) Charge pumped per cycle from left 
to right versus the Fermi wave number for a static potential 
separated from an oscillating potential by five sites, with o = 
1, b = 0.25, and uj — tt/10; the system has 252 sites. The 
numerical and analytical results with and without the bound 
state contribution are shown by dash-dot (red), starred (blue) 
and dashed (black) lines, respectively. 



FIG. 4: (Color online) Charge pumped per cycle from left 
to right versus the oscillation frequency for a static potential 
separated from the oscillating potential by five sites, with 
a = 1.5, b — 0.25, and &f = 7i"/4; the system has 308 sites. 
The numerical and analytical results are shown by dash-dot 
(red) and starred (blue) lines, respectively. 



dot and starred lines show the numerical and analytical 
results obtained from Eq. (fTTj) respectively. A noticeable 
change is seen to occur when to crosses a value of about 
0.59; this is because the first side band E-\ = /i — u> 
goes below the continuum band and stops contributing 
to the pumped charge at that point. A less noticeable 
change occurs for a similar reason when uj crosses a value 
of about 3.41, where the first side band E\ = /i + uj 
goes above the continuum band. We also note that the 
pumped charge goes to zero at u> = 4, as we had com- 
mented after Eq. pip . 

Finally, Fig. 5 shows the charge pumped per cycle 
as a function of the oscillation amplitude b for a = 1.5, 
uj = 7r/10, and kp — -k/A. In this case, we cannot use 
the analytical expression given in Eq. (jlip since b is not 
small in general and there is a substantial contribution 
from higher side bands with \p\ > 2. It is interesting to 
note that the pumped charge vanishes and changes sign 
at certain values of b. 



IV. DISCUSSION 

In this paper, we have studied charge pumping by a sin- 
gle oscillating potential when the spatial left-right sym- 
metry is broken due to the presence of a static potential 
at another point; this model is similar to the one used to 
describe some recent experiments [38] • We have used Flo- 
quet scattering theory to compute the pumped charge to 
second order in the oscillation amplitude. We have shown 
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FIG. 5: (Color online) Charge pumped per cycle from left 
to right versus the oscillation amplitude for a static potential 
separated from the oscillating potential by five sites, with 
a — 1.5, uj — tt/10, and /cf = tt/4; the system has 252 sites. 



that if the oscillation frequency is larger than a threshold 
value, a bound state can contribute to charge pumping; 
this possibility does not seem to have been studied ear- 
lier. 

For small amplitudes, we find that the results obtained 
numerically using the equation of motion method gener- 
ally agree well with the analytical results. However, some 
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numerical problems arise when the contribution from a 
bound state becomes important. These problems have 
been observed earlier in Ref. [5l| (see also Ref. [EH); 
they are due to a difficulty in maintaining the occupa- 
tion of the bound state at the correct equilibrium value. 
The simple model we have used for the numerical cal- 
culations has no interactions between electrons and no 
phonons which can lead to inelastic scattering processes 
and thereby maintain the occupation of the bound state 
at a value dictated by the Fermi function. Some ways of 
addressing the problem of bound states have been dis- 
cussed in Refs. [HlllSa ]. 

It would be interesting to generalize our analysis to the 
case of an arbitrary time-dependent potential applied at 
one point where the potential may contain a very large 
number of oscillation frequencies with arbitrary ampli- 
tudes. In particular, one can consider the case of a noisy 



potential and study whether that can, on the average, 
pump charge in one direction when the left-right symme- 
try is broken by a static potential. The brief discussion in 
Sec. II of a potential consisting of a few oscillating terms 
with low frequencies and small amplitudes suggests that 
a weak noise may indeed be able to pump charge on the 
average, but a detailed investigation of this may be use- 
ful. 



Acknowledgments 

A. A. thanks CSIR, India for financial support. D.S. 
thanks A. Dhar for many stimulating discussions. We 
thank DST, India for financial support under the projects 
SR/FST/PSI-022/2000 and SR/S2/CMP-27/2006. 



[1] D. J. Thouless, Phys. Rev. B 27, 6083 (1983). 

[2] Q. Niu, Phys. Rev. Lett. 64, 1812 (1990). 

[3] M. Biittiker, H. Thomas, and A. Pretre, Z. Phys. B: Con- 
dens. Matter 94, 133 (1994). 

[4] P. W. Brouwer, Phys. Rev. B 58, R10135 (1998). 

[5] P. W. Brouwer, Phys. Rev. B 63, 121303(R) (2001). 

[6] I. L. Aleiner and A. V. Andreev, Phys. Rev. Lett. 81, 
1286 (1998). 

[7] C. S. Tang and C. S. Chu, Solid State Commun. 120, 
353 (2001). 

[8] S. W. Chung, C. S. Tang, C. S. Chu, and C. Y. Chang, 

Phys. Rev. B 70, 085315 (2004). 
[9] J. E. Avron, A. Elgart, G. M. Graf, and L. Sadun, Phys. 

Rev. B 62, R10618 (2000); Phys. Rev. Lett. 87, 236601 

(2001) ; J. Stat. Phys. 116, 425 (2004). 

[10] O. Entin-Wohlman, A. Aharony, and Y. Levinson, Phys. 

Rev. B 65, 195411 (2002). 
[11] O. Entin-Wohlman and A. Aharony, Phys. Rev. B 66, 

035329 (2002). 

[12] N. Y. Hwang, S. C. Kim, P. S. Park, and S.-R. Eric Yang, 

larXiv:0706.0947l 
[13] Y. M. Galperin, O. Entin-Wohlman, and Y. Levinson, 

Phys. Rev. B 63, 153309 (2001). 
[14] A. Aharony and O. Entin-Wohlman, Phys. Rev. B 65, 

241401 (R) (2002). 
[15] V. Kashcheyevs, A. Aharony, and O. Entin-Wohlman, 

Phys. Rev. B 69, 195301 (2004), and Eur. Phys. J. B 39^ 

385 (2004). 

[16] M. Moskalets and M. Biittiker, Phys. Rev. B 66, 205320 

(2002) , and Phys. Rev. B 68, 075303 (2003). 
[17] S. W. Kim, Phys. Rev. B 66, 235304 (2002). 

[18] B. Wang, J. Wang, and H. Guo, Phys. Rev. B 65, 073306 
(2002), and Phys. Rev. B 68, 155326 (2003). 

[19] L. Arrachea, Phys. Rev. B 72, 125349 (2005), and Phys. 
Rev. B 72, 121306(R) (2005). 

[20] L. E. F. Foa Torres, Phys. Rev. B 72, 245339 (2005). 

[21] L. Arrachea, C. Naon, and M. Salvay, Phys. Rev. B 76, 
165401 (2007). 

[22] A. Banerjee, S. Das, and S. Rao, 

|arXiv:cond mat /0307324| 
[23] S. Banerjee, A. Mukherjee, S. Rao, and A. Saha, Phys. 



Rev. B 75, 153407 (2007). 
[24] G. R. Aizin, G. Gumbs, and M. Pepper, Phys. Rev. B 58, 

10589 (1998); G. Gumbs, G. R. Aizin, and M. Pepper, 

Phys. Rev. B 60, R13954 (1999). 
[25] K. Flensberg, Q. Niu, and M. Pustilnik, Phys. Rev. B 

60, R16291 (1999). 
[26] P. A. Maksym, Phys. Rev. B 61, 4727 (2000). 
[27] A. M. Robinson and C. H. W. Barnes, Phys. Rev. B 63, 

165418 (2001). 

[28] M. M. Mahmoodian, L. S. Braginsky, and M. V. Entin, 

Phys. Rev. B 74, 125317 (2006). 
[29] M. Strass, P. Hanggi, and S. Kohler, Phys. Rev. Lett. 95, 

130601 (2005). 

[30] S. Kohler, J. Lehmann, and P. Hanggi, Phys. Rep. 406, 
379 (2005). 

[31] A. Agarwal and D. Sen, J. Phys. Condens. Matter 19, 
046205 (2007). 

[32] M. Switkes, C. M. Markus, K. Campman, and A. C. 

Gossard, Science 283, 1905 (1999). 
[33] V. I. Talyanskii, J. M. Shilton, M. Pepper, C. G. Smith, 

C. J. B. Ford, E. H. Linfield, D. A. Ritchie, and G. A. C. 
Jones, Phys. Rev. B 56, 15180 (1997). 

[34] J. Cunningham, V. I. Talyanskii, J. M. Shilton, M. Pep- 
per, M. Y. Simmons, and D. A. Ritchie, Phys. Rev. B 
60, 4850 (1999); J. Cunningham, V. I. Talyanskii, J. M. 
Shilton, M. Pepper, A. Kristensen, and P. E. Lindelof, 
Phys. Rev. B 62, 1564 (2000). 

[35] V. I. Talyanskii, D. S. Novikov, B. D. Simons, and L. S. 
Levitov, Phys. Rev. Lett. 87, 276802 (2001). 

[36] P. J. Leek, M. R. Buitelaar, V. I. Talyanskii, C. G. Smith, 

D. Anderson, G. A. C. Jones, J. Wei, and D. H. Cobden, 
Phys. Rev. Lett. 95, 256802 (2005). 

[37] L. DiCarlo, C. M. Marcus, and J. S. Harris, Jr., Phys. 

Rev. Lett. 91, 246804 (2003); M. G. Vavilov, L. DiCarlo, 

and C. M. Marcus, Phys. Rev. B 71, 241309(R) (2005). 
[38] B. Kaestner, V. Kashcheyevs, S. Amakawa, L. Li, M. 

D. Blumenthal, T. J. B. M. Janssen, G. Hein, K. 

Pierz, T. Weimann, U. Siegner, and H. W. Schumacher, 

larXiv:0707.0993l 
[39] S. Das and S. Rao, Phys. Rev. B 71. 165333 (2005). 
[40] C. de C. Chamon, D. E. Freed, and X. G. Wen, Phys. 



8 



Rev. B 51, 2363 (1995) and Phys. Rev. B 53, 4033 (1996). 
[41] C. dc C. Chamon, D. E. Freed, S. A. Kivelson, S. L. 

Sondhi, and X. G. Wen, Phys. Rev. B 55, 2331 (1997). 
[42] P. Sharma and C. Chamon, Phys. Rev. Lett. 87, 096401 

(2001), and Phys. Rev. B 68, 035321 (2003). 
[43] D. E. Feldman and Y. Gefen, Phys. Rev. B 67, 115337 

(2003). 

[44] D. Makogon, V. Juricic, and C. Morais Smith, Phys. 
Rev. B 74, 165334 (2006), and Phys. Rev. B 75, 045345 
(2007). 

[45] D. Schmeltzer, Phys. Rev. B 63, 125332 (2001). 

[46] D. E. Feldman, S. Scheidl, and V. M. Vinokur, Phys. 



Rev. Lett. 94, 186809 (2005). 
[47] F. Cheng and G. Zhou, Phys. Rev. B 73, 125335 (2006). 
[48] B. Braunecker, D. E. Feldman, and F. Li, Phys. Rev. B 

76, 085119 (2007). 
[49] T. L. Schmidt and A. Komnik, Phys. Rev. B 76, 035113 

(2007). 

[50] A. Agarwal and D. Sen, Phys. Rev. B 76, 035308 (2007). 
[51] A. Dhar and D. Sen, Phys. Rev. B 73, 085119 (2006). 
[52] A. Agarwal and D. Sen, Phys. Rev. B 73, 045332 (2006). 
[53] G. Stefanucci, Phys. Rev. B 75, 195115 (2007). 



